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ABSTRACT 

Substellar Objects in Nearby Young Clusters - SONYC - is a survey program to investigate the 
frequency and properties of substellar objects in nearby star-forming regions. We present new spec¬ 
troscopic follow-up of candidate members in Chamaeleon-I 2Myr, 160 pc) and Lupus 3 (^IMyr, 

200pc), identified in our earlier works. We obtained 34 new spectra (1.5 - 2.4 /rm, R^600), and iden¬ 
tified two probable members in each of the two regions. These include a new probable brown dwarf 
in Lupus 3 (NIR spectral type M7.5 and Teg = 2800 K), and an L3 {Tes = 2200 K) brown dwarf in 
Cha-I, with the mass below the deuterium-burning limit. Spectroscopic follow-up of our photometric 
and proper motion candidates in Lupus 3 is almost complete (> 90%), and we conclude that there are 
very few new substellar objects left to be found in this region, down to 0.01 - 0.02Mq and Ay < 5. 

The low-mass portion of the mass function in the two clusters can be expressed in the power-law form 
dN/dM (X M““, with a ^ 0.7, in agreement with surveys in other regions. In Lupus 3 we observe a 
possible flattening of the power-law IMF in the substellar regime: this region seems to produce fewer 
brown dwarfs relative to other clusters. The IMF in Cha-I shows a monotonic behavior across the 
deuterium-burning limit, consistent with the same power law extending down to 4 - 9 Jupiter masses. 

We estimate that objects below the deuterium-burning limit contribute of the order 5 — 15% to the 
total number of Cha-I members. 

Subject headings: stars: formation, low-mass, brown dwarfs, mass function 


1. INTRODUCTION 

SONYC - short for Substellar Objects in Nearby Young 
Clusters - is a comprehensive project aiming to provide a 
complete, unbiased census of substellar population down 
to a few Jupiter masses in young star forming regions. 
The survey is based on extremely deep optical- and near- 
infrared wide-held imaging, combined with the Two Mi¬ 
cron All Sky Survey (2MASS) and Spitzer photometry 
catalogs, which are correlated to create catalogs of sub¬ 
stellar candidates and used to identify targets for exten¬ 
sive spectroscopi c follow-up. To furt her facilitate candi¬ 
date selection, in Muzic et al. (20141 we also performed 
a proper motion analysis of the Lupus 3 star forming 
region. While the candidate selection based on optical 
and near-infrared photometry helps to avoid biases in¬ 
troduced by the mid-infrared selection (only objects with 
disks), or methane-imaging (only T-dwarfs), it comes 
with a cost of relatively large candidate lists affected 
by significant background contamination. Spectroscopic 
follow-up is therefore a necessary prerequisite to reliably 
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characterize the low mass population in young star form¬ 
ing regions. 

Thanks to our SONYC survey and the efforts of other 
groups, the substellar Initial Mass Function (IMF) is now 
well characterized down to ^0.01 Mq. The ratio of the 
number of stars with respect to b rown dwarfs lies be¬ 
tween 2 and 6 (Scholz et al. 20131, and the power-law 
slope of the mass function is a ~0.6 (dN/dMocM~°' 


see 


Scholz et al. 2012a Offner et al. |26l4| for a sum- 


mary^ It is clear by now that the mass function in 
young clusters extends below the deuterium-burning (D- 
burning) limit at ~ 13Mj„p, as a handful of early L- 


NGC1333 (IScholz et al.|2012a), p Oph I 

Alves de Oliveira 

et al. 20121), A Uri (Bayo et al. 

2011 

1 , a Ori (IZapa- 

tero Osorio et al. 

200Dr2013), Orion 

Ingraham et al. 

20I4|, and Cha-I 

ILuhman et al. 

2008 

. However, the 


^ lOMjup, leaving the mass function in the planetary- 
mass regime still poorly constrained. The surveys that 
extend to lower masses are mostly based on photome¬ 
try, and still await spectroscopic follow-up. The SONYC 
study in NGC 1333 (Scholz et al. 2012b) is unique in 
this sense, as we obtained spectra of more than 85% of 
all photometric very-low-mass (VLM) candidates down 
to < O.OO5M0. We find that the free-floating objects 
with planetary masses are rare, 20-50 times less numer¬ 
ous than stars (for planetary mass objects in range 0.006 
- O.O15M0, and stars below IM 0 ). The mass spectrum 
in NGC 1333 below 0.6 M 0 is well described by a power- 
law form dN/dM oc M““, with a « 0.6, and possibly 
requires l ower values of a in the plan etary-mass domain. 
In cr Ori, Pena Ramirez et al. (2012) find the same slope 
for the mass interval 0.25-0.004 M 0 , and also suggest a 
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possibly lower a below 0.004Mq. While this result still 
awaits spectroscopic confirmation, three T-dwarf candi¬ 
dates tow ards a Ori have been classif ied as likely non¬ 
members (Pena Ramirez et al. 20151, suggesting that 
the mass function in this cluster might have a cut-off 
at ^4 Mjup. In the central part of Upper Scorpius , 
spectroscopy of L-type candidates (iLodie u et al.||2011|) 
and pr oper motions of T-type candidates (jbodieu et al. 
2013b I favor a turn-down of the mass function below 10 
- 4Mj„p, depending on the age assigned to the cluster, 
and the models used. In another photometric survey 
probing the planetary-ma ss range in th e same region, 
but over much larger field, 
the mass function (in the 


Lodieu (2013) concludes that 
og-norrnal form) is likely de¬ 
creasing in the planetary-mass regime, although the flat 
function cannot be discarded at this point. A proper- 
motion and/or spectroscopic analysis of the candidate 
members is needed before drawing a firmer conclusion. 

The surveys in star forming regions thus far are possi¬ 
bly hi_conhi£twith_the claims of the microlensing study 
by|Sumi et al.| (|2011|), who suggest that the unbound or 
wide-orbit Jupiter-mass objects are almost twice as com¬ 
mon as main sequence stars. One possible mechanism for 
the formation of objects of this kind might be ejection of 
BDs, and/or star forming clumps through close encoun¬ 
ters in the early stages of cluster formation. However, the 
typical ejection velocities produced by simulations (e.g. 


Bate 

120091120121 IBasu & Vorobyov||2012l 

Reipurth et al.j 

2010 

W-H] 

) of up to a few kms“^ might result 
Ds distributed at the outskirts of tt 

n a population 
le clusters, but 


not very far from it. Certainly some objects at the tail of 
the velocity distribution might escape even further, but 
this cannot be a substantial number. For example, at the 

rlB — 

BI 


end of the cluster simulation in|Bate (2009), only about 
17% of the low-mass stars and BDs are at the distances 
beyond the 80%-mass radius of the cluster. Therefore, 
if the microlensing result holds, the objects they probe 
should have masses below ~ 0.005 Mq, and their over¬ 
abundances would favor a different formation scenario 
from the one forming the non-D-burning brown dwarfs 
in young clusters. 

Within SONYC, we have investigated the IMF in the 
two yo ung star forming reg ions, Chamaeleo n-I (hereafter 


Cha-I ; Muzic et al. 2011), and Lupus 3 (Muzic et al. 


2014). In this paper we present the second spectroscopic 


follow-up of the SONYC candidates in these two regions. 
Based on this more complete set of follow-up spectra, 
we can better constrain the properties of the IMF at 
the low-mass end. In Cha-I specifically, we target the 
candidates expected to have masses below the D-burning 
limit, in order to derive more definitive conclusions that 
contribute to the ongoing discussion about the IMF in 
this mass regime. 

In Section we give a summary of our previous work 
in the two regions, describing in detail the candidate lists 
and the first spectroscopic follow-up. In Section we 
present the new spectroscopic follow-up with Sofl/NTT. 
The analysis of the spectra is presented in Section]^ The 
results are discussed in Section [^Finally, we summarize 
the main conclusions in Section 

2. CANDIDATE LISTS AND SPECTROSCOPIC FOLLOW-UP 

In this section we outline p revious SONY C campaigns 
in Cha-I ( Muzic et al.||20lT ), and Lupus 3 ( [Muzic et al. 


TABLE 1 

Overview of the various samples in Cha-I and Lupus 3 
USED IN THIS paper 


Sample 

No.*" 

Chamaeleon-I 


Candidates selected from Iz diagram 
with spectra from VIMOS 
with spectra from SofI 
confirmed VLM objects by VIMOS spectra 
confirmed VLM objects by SofT spectra 
confirmed VLM objects by other groups'^ 
rejected as VLM objects by our spectra 
rejected as VLM objects by other groups^ 

142/106 

18/18 

15/15 

7/7 

2/2 

9/9 

24/24 

4/4 

Lupus 3 


Candidates selected from ij diagram^ 
with spectra from VIMOS 
with spectra from SofI 

confirmed young VLM objects by VIMOS spectra 
confirmed young VLM objects by SofI spectra 
confirmed young VLM objects by other groups^ 

372/337 

123/112 

19/19 

7/7 

2/2 

1/1 

Candidates from “IJ-pm” sample^ 
with spectra from VIMOS 
with spectra from SofI 

confirmed young VLM objects by VIMOS spectra 
confirmed young VLM objects by SofI spectra 
confirmed young VLM objects by other groups 

58/53 

32/31 

19/19 

7/7 

2/2 

0/0 


^ total number/above the completeness limit 

^ in addition to the members confirmed/rejected by our spectra 
^ all photometric candidates, including the “IJ-pm” sample 
^ ij and proper motion candidates 


2014), as the analysis presented in this work is a con- 
tinuation of our past efforts. We give a summary of the 
candidate lists for the two regions, as well as a detailed 
description of the second follow-up with Sofl/NTT. In 
Table we provide an overview of the candidate samples 
used fm' the analysis presented in Section [^ For each en¬ 
try, we give the total number of objects, and the number 
of objects that are brighter than the completeness limit 
in each of the two surveys. The latter is important since 
we can derive firm statistical conclusions on the shape 
of the IMF only in the region where our surveys can be 
considered complete. 


2.1. Chamaeleon-I 

While the Cha-I star forming region spreads over sev¬ 
eral square degrees on the sky, more than 80% of the 
known membe rs lie within a ~ 1 deg^ central stripe (see 
Luhman 2007). The SONYC photometric survey pro- 
vided deep images of a part of this area, covering in total 
~0.25deg^. Our fields contain 38% of the Cha-I pop¬ 
ulation accor ding to the mo st complete census to date 
presented in Luhman (2007). Our photometric analy¬ 
sis was based on optical and near-infrared catalogs in I, 
z, J, and Ks band^ A list of candidate members of 
Cha-I was selected from the (/ — z, I) color-magnitude 
diagram (CMD), representing a natural extension of the 
previously known members towards fainter magnitudes. 


® Throughout this paper, I stands for Cousins I, and i for DENIS 

1- band passband. NIR photometry is in the 2MASS system. The 

2 - band photometry was not absolutely calibrated, and is only used 
to s eparate the red fro m the blue sources in the {I — z) color space 
(seejMuzic et al.|2011|for details). 
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Fig. 1.— Left: Color-magnitude diagram showing candidate members of Cha-I (open circles). The filled circles represent those also found 
in our NIR catalog. The sources observed in two spectroscopic campaigns are shown as black squares (VIMOS/VL T), and blue tria ngles 
(Sofl/NTT). The confirmed members of Cha-I are marked with red diamonds, and those classified as non-members in|Luhman|l|2004|| with 
crosses. The horizontal lines in the left panel mark the upper limit for the Soil follow-up (at I = 21), and the 7-band completeness limit 
(7 = 23). Right: Color-magnitude diagram showing candidate members in Lupus 3 (dots). Proper motion candidates are additionally 
marked with open circles. The remaining symbols are the same as in the left panel, la photometric uncertainties are shown on the left-hand 
side of the plot. The expected object masses at Ay = 0 according to the BT-Settl models are shown on the right-hand side of each plot, 
for an age of 2 Myrs in Cha-I, and 1 Myr in Lupus 3. 


This list contains 142 candidate members, with /-band 
magnitudes between 17.5 and 24. For Cha-I m embers 
within this brightness range, BT-Settl models (Allard 
et al. 2011) predict masses roughly between 0.002 and 
O.OOMq, for Ay < and an age of 1-2 Myr. The com¬ 
pleteness limit of the /—band catalog is 23 mag, which 
is equivalent to ^ 0.003 — 0.007 Mq for Ay < 5 at 1 
Myr, or ~ 0.004 — O.OOOMq at the age of 2 Myrs (me¬ 
dian age of Cha-I members when compared to theoretical 
isochrones; Luhman 2007), and the distance of 160 pc. 
The completeness limit of the near-infrared catalogs is 
J = 18.3 and Ks = 16.7, which corresponds to masses 
- 0.004 - 0.007 Mq for Ay < 5 
The first spectroscopic follow-up of the objects in the 
candidate list was performed using the multi-object spec¬ 
trograph VIMOS at the VLT. The follow-up included 
18 objects from the candidate list, most of which have 
/ < 21. We confirmed 7 objects as VLM members of 
Cha-I. In the left panel of the Figure ^ we show the 
CMD for candidate members. For the full CMD we re- 


Muench (12008), while crosses mark non-members from 


fer the reader to Figure 3 in Muzic et al. (2011). Two 
horizontal lines in the left panel mark the approximate 
limit of the VIMOS follow-up (/ = 21), and the I- 
band completeness limit (/ = 23). If indeed members of 
Cha-I, the objects between these two lines should have 
masses between 0.003 and 0.015 Mq at IMyr or between 
0.004 and 0.02 Mq at 2 Myr (for Ay < 5). Red di¬ 
amonds mark spectroscop ically con f irmed members of 
Cha-I, from SONYC and Luhman (2007); Luhman & 


^ According to the extinction maps of ICambrfevI jl999ll. only 

abou^J_%_^mtieUIiaUa£eaUiusA^_>5jWhile_^iieeiArej/u£us_3 

a£ea_amears_to_have^j^<5^_We_^heretoreuse_^his_^san_u££er 

value ot extinction when deriving mass limits irom the photometry. 


Luhman (|2UU4). 

For the second follow-up with the Sofl/NTT, we tar¬ 
geted the sources with expected masses in the planetary- 
mass regime detected in our NIR images (the filled dots 
in the left panel of Figurej^ located between the two hor¬ 
izontal lines), and not previously observed with VIMOS 
(squares). The 15 SofI follow-up targets are marked with 
blue triangles, and listed in Table 

2.2. Lupus 3 

The SONYC photometric survey in Lupus 3 covers 
1.4 deg^ of the main cluster core, surrounding the two 
brightest members H R599 9/6000. Previous s urveys by 
Comeron et al. (2009) and Merfn et al. (2008) cover ap- 
proximately this same area, plus the lower-density part 
of the complex to the North-East of the main core. This 
part, however, contains only about 10% of all the can¬ 
didate members identified in the two aforementioned 
works. The original list of photometric candidates in 
Lupus 3 presented in Muzic et al. (2014) was selected 
from the {i — J, i) CAlU, and contained 409 candidates. 
Although in our catalogs we removed the sources too 
close to the edges of the chips, and overly elongated or 
blended sources (based on SExtractor output), some ex¬ 
tended sources (galaxies with almost round profiles, un¬ 
recognized doubles) and a few artifacts (source sitting 
right on a bad column, or a spike of a nearby bright 
object) still made it to our candidate catalog. Total of 
37 contaminants of this type were found by visual in¬ 
spection of the images, and have been removed from the 
candidate list. The final candidate list then contains 372 
sources. This list can be further narrowed-down to 58, 


taking into account the proper motions (see Muzic et al. 
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2014 

this 


for details). As in the previous paper, we refer to 
ist as the “IJ-pm” candidate list. The photometric 


completeness limit of our catalog is i=20.3, equivalent to 
0.008 — 0.02 Mq for Ay = 0 — 5 at a distance of 200 
pc and age of 1 Myr, according to the BT-Settl models. 
The magnitude range of the candidate list / = 15.3 — 21 
is eq uivalent to masses 0. 2 - O.OOTMq for Ay < 5. 

In Muzic et al. (2014), we presented the first spec- 
troscopic follow-up using VIMOS/VLT, where we took 
spectra of 123 photometric candidates, out of which 32 
are from the “IJ-pm” list. Spectra confirmed 7 of these 
objects as young members of Lupus 3, all belonging to 
the “IJ-pm” list. For the second follow-up, the SofI spec¬ 
troscopy presented here, the targets were selected from 
the “IJ-pm” list not covered by the VIMOS follow-up, 
and above the completeness limit. The second follow-up 
includes 19 candidates, listed in Table 

The right panel of the Figure shows the CMD con¬ 
taining 372 candidate members of Lupus 3 (dots). The 
proper motion candidates (“IJ-pm”) are additionally 
marked with open circles. Spectroscopic targets from 
the two follow-ups (VIMOS and SofI) are marked with 
black squares and blue triangles, respectively. 

2.3. Spectroscopy with Sofl/NTT 

The spectra were ta ken using SofI (Son of ISAAC; 
Moorwood et al.||1998|) at the ESO’s New Technology 
'i'elescope (IN'IT’) during four nights in April 2014, un¬ 
der the program number 093.C-0050(A). SofI is equiped 
with a Hawaii HgCdTe 1024 x 1024 detector, with the 
pixel size projected onto the sky of 0.288". We used 
the low resolution red grism (GR Grism Red) covering 
the wavelength range 1.5 - 2.5 /rm, and delivering spec¬ 
tral resolution R~ 600 when combined with the 1" slit 
(31 spectra), and 300 with the 2" slit (3 spectra). The 
seeing during the SofI observations was varying between 
0.7 and 1.5" (values from the Differential Image Motion 
Monitor; DIMM), and the transparency was clear, ex¬ 
cept for the first few hours of the night 2014-04-17, when 
some clouds have been present. The wider slit was used 
for the patricularly faint sources {K ^ 16), during less 
favorable seeing conditions. The typical nod throw along 
the slit was 60". The complete SofI target list, along with 
the information about the date of observation, on-source 
integration time, and slit used, is presented in Table 

Data reduction was performed by combining the SOFI 
pipeline recipes, our IDL routines, and IRAF task apall. 
The data reduction steps include cross-talk, flat field, and 
distortion corrections. Pairs of frames at different nod¬ 
ding positions were subtracted one from another, shifted, 
and combined into a final frame. The spectra were ex¬ 
tracted using the IRAF task apall, and wavelength cali¬ 
brated with the help of the Xenon lamp arcs. The wave¬ 
length solution obtained by the pipeline has an RMS 
error of ~ 20A. The spectra of telluric standard stars 
were reduced in the same way as science frames. Telluric 
standards spectral types range between B5V and B9.5V, 
and show several prominent hydrogen lines in absorption, 
which were removed from the spectra by interpolation, 
prior to division with the black body spectrum at an ap¬ 
propriate effective temperature. This yields the response 
function of the spectrograph, convolved with the telluric 
spectrum. Finally, we divided the science spectra by the 
corresponding response function. The error of the re¬ 


sponse function is the combination of the errors of the 
measured telluric spectra and the errors in the effective 
temperature assigned to certain spectral types. The lat¬ 
ter is caused by assumption that the instrinsic spectrum 
of the calibrating star is well described by a black body 
of a certain temperature. The error of the measured tel¬ 
luric spectra is supplied by apall, through the extras 
keyword. It is of the order, or below 1% across most of 
the spectrum, except in the region between H and K- 
bands where the atmosphere absorbs most of the light. 
For the B-type tell uric standards, we adop ted the tem¬ 
perature scale from Boehm-Vitense (1981), whose ±5% 
error on the effective temperature translates to approxi¬ 
mately ±1% error in the response curve. 

3. DATA ANALYSIS 

The analysis of the spectra presented in this section 
consists of: (1) visual inspection of the spectra in search 
for the features revealing the youth, (2) fitting of the 
model spectra to the observed ones to determine the ef¬ 
fective temperature and extinction, and (3) determina¬ 
tion of spectral types. 

3.1. Membership assessment 

The initial assessment of the spectra is done by vi¬ 
sual inspection. Spectra of young objects later than 
' M5 show a characteristic triangu lar peak in the H- 


band, caused by w ater absorption (Lucas et al. 2001 
Gushing et al.||2005 ). The depth of this feature strongly 
depends on ettective temperature, while the shape is af¬ 
fected by the gravity. The peak appears triangular for 
young, low- to intermediate-gravitjj^objects, and round 
at higher gravities. The K-band spectra of late type ob¬ 
jects show prominent GO absorption bands. GO bands 
can facilitate spectral classification, since K and M-type 
giants show (deeper GO absorption than dwarfs (see e.g. 
Rayner et al.||2069 ). We also look for the emission lines, 
which may indicate accretion. Only three objects from 
the SofI spectroscopic sample (1 in Gha-I, 2 in Lupus 3) 
clearly show a triangular H-band peak. They are shown 
in Table and Figure (black), along with the cor¬ 
responding best fit modms (red; see the next section). 
SONYG-Lup3-29 is a new, probably sub-stellar member 
of Lupus 3. The membership is further confirmed by the 
mid-infrared excess, as demonstrated in Figure!^ 

In Figure we also show the spectrum of ISO-Ghal- 
192, whose steep spectrum with prominent H 2 emission 
lines indicates a YSO nature Indeed, this objec t is a 
known class I protostar and an FU Ori candidate ( Gra- 


majo et al. 12014 

Gambresv et al.|1998 

Persi et al.|1999D, 

associated with a GO outflow (|Mattila et al. 

TUS^rP^sil 


type ot lV13.5dVI6.5 from the near infrared low-resolution 
spectra. Two estimates for the ef fective temperatur e are 
found in t he literature, 3600K b y Persi et al.] ( 2007[ ), and 


5000 K by Gramajo et al. (2014). i'he typical class-i SED 


of ISO-GhaI-192 is shown in Figure [3 
The four probable members and their properties de¬ 
rived in the next sections are listed in Table I The re- 


< 300 Myr according to 


Allers & Liu 


(20131. 


The clearly detected H 2 lines are 1-0 y(2) at 2.034^m, 1-0 

S(l) at 2.122/im, 1-0 S(3) at l.gST/rm, 1-0 S(0) at 2.222/rm, 2-1 

S(l) at 2.247/rm, 1-0 Q(l) at 2.404/rm, 1-0 Q(3) at 2.421/rm, and 

1-0 Q(4) at 2.435Atm. 
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TABLE 2 

Photometric candidates included in the spectroscopic follow-up with SofI 


# 

o(J2000) 

(5(J2000) 

Date 

Exp. time Slit 


J 

K 

comments 

Chamaeleon-I 


1 

11 

05 

21.36 

-77 

46 

33.0 

2014-04-19 

14 X 180 s 

1" 

21.92 

17.37 

14.65 


2 

11 

06 

18.62 

-77 

35 

17.4 

2014-04-18 

12 X 180 s 

1" 

22.43 

16.70 

13.22 


3 

11 

06 

28.04 

-77 

35 

54.7 

2014-04-13 

14 X 60s 

1" 

21.68 

16.50 

13.31 


4 

11 

08 

20.20 

-77 

45 

48.6 

2014-04-13 

10 X 200 s 

1" 

21.50 

18.09 

15.86 


5 

11 

08 

23.59 

-77 

31 

29.5 

2014-04-14 

52 X 100 s 

2" 

22.26 

18.65 

16.49 


6 

11 

08 

30.31 

-77 

31 

38.6 

2014-04-14 

52 X 100 s 

2" 

22.50 

17.92 

16.05 

confirmed member 

7 

11 

08 

30.77 

-77 

45 

50.5 

2014-04-13 

10 X 200 s 

1" 

21.34 

17.77 

15.59 


8 

11 

08 

33.92 

-77 

46 

33.4 

2014-04-13 

14 X 220 s 

1" 

21.49 

18.54 

16.68 


9 

11 

08 

50.38 

-77 

45 

53.0 

2014-04-13 

14 X 220 s 

1" 

21.26 

17.90 

16.01 


10 

11 

09 

11.62 

-77 

19 

26.1 

2014-04-18 

10 X 180 s 

1" 

21.74 

16.56 

13.20 


11 

11 

09 

28.57 

-76 

33 

28.1 

2014-04-14 

18 X 90s 

1" 

21.95 

16.74 

11.93 

confirmed member 

12 

11 

09 

32.26 

-76 

34 

39.3 

2014-04-19 

20 X 90s 

2" 

22.76 

18.52 

15.74 


13 

11 

09 

35.44 

-77 

17 

27.1 

2014-04-18 

8 X 180 s 

1" 

21.68 

16.73 

13.60 


14 

11 

09 

37.38 

-76 

41 

05.7 

2014-04-17 

16 X 250 s 

1" 

22.71 

17.92 

15.18 


15 

11 

09 

51.82 

-76 

42 

18.2 

2014-04-19 

12 X 180 s 

1" 

22.97 

17.93 

14.57 


Lupus 3 

16 

16 

08 

03.96 

-39 

10 

28.9 

2014-04-18 

10 X 150 s 

1" 

19.39 

16.28 

14.21 


17 

16 

08 

25.87 

-39 

10 

08.0 

2014-04-17 

8 X 120 s 

1" 

18.13 

15.26 

13.03 


18 

16 

08 

30.04 

-39 

10 

01.4 

2014-04-17 

8 X 120 s 

1" 

15.51 

12.69 

10.82 


19 

16 

08 

33.04 

-38 

52 

22.7 

2014-04-18 

6 X 180 s 

1" 

15.34 

12.98 

11.74 

confirmed member 

20 

16 

09 

42.66 

-39 

03 

20.7 

2014-04-18 

10 X 180 s 

1" 

18.68 

15.54 

13.29 


21 

16 

09 

43.74 

-39 

07 

11.9 

2014-04-17 

12 X 120 s 

1" 

18.96 

15.74 

13.23 


22 

16 

09 

47.20 

-39 

08 

42.4 

2014-04-19 

6 X 180 s 

1" 

18.92 

15.91 

13.82 


23 

16 

09 

52.68 

-39 

06 

40.4 

2014-04-17 

8 X 120 s 

1" 

17.13 

14.45 

12.52 


24 

16 

10 

01.33 

-39 

06 

45.1 

2014-04-17 

8 X 60s 

1" 

15.38 

12.32 

10.52 

confirmed member 

25 

16 

10 

05.97 

-39 

08 

06.2 

2014-04-18 

10 X 180 s 

1" 

18.40 

15.54 

13.29 


26 

16 

10 

06.01 

-39 

04 

23.7 

2014-04-17 

6 X 60s 

1" 

16.90 

12.30 

9.59 

giant? 

27 

16 

10 

11.71 

-39 

09 

19.8 

2014-04-17 

8 X 180 s 

1" 

17.81 

14.84 

12.96 


28 

16 

10 

13.74 

-39 

02 

34.8 

2014-04-17 

10 X 120 s 

1" 

17.21 

13.49 

10.86 

giant? 

29 

16 

10 

14.59 

-39 

02 

31.4 

2014-04-17 

10 X 120 s 

1" 

19.41 

15.89 

13.58 


30 

16 

10 

17.86 

-39 

03 

47.1 

2014-04-17 

6 X 120s 

1" 

16.26 

13.73 

11.93 


31 

16 

10 

30.82 

-39 

04 

04.4 

2014-04-18 

6 X 180 s 

1" 

17.29 

14.66 

12.81 


32 

16 

10 

33.59 

-39 

08 

21.4 

2014-04-17 

6 X 180 s 

1" 

17.56 

14.84 

12.54 


33 

16 

11 

46.57 

-39 

06 

04.7 

2014-04-17 

8 X 90s 

1" 

15.39 

13.19 

11.23 

member? 

34 

16 

12 

04.95 

-39 

06 

25.6 

2014-04-17 

6 X 120 s 

1" 

15.40 

13.12 

11.51 

giant or member 


^ Cousins I for sources in Chamaeleon, and DENIS i for those in Lupus 


maining objects show flat spectra with few features, and 
cannot be easily classified at the resolution and signal-to- 
noise of our spectra. They are shown in Figures 15 and 
for Lupus 3, and Cha-I, respectively, and listed in Table^ 
(all objects in this table except those with the comment 
“confirmed member”). The fact that their spectra do 
not show the water absorption in the H-band indicative 
of cool spectral types means that they must have spec¬ 
tral type earlier than coM5, i.e. they are certainly not 
substellar. The lack of CO bands at 2.3/im in most of 
these spectra is another signpost of warm atmospheres. 
Only the top-most four spectra of the right panel in Fig¬ 
ure clearly show CO absorption bands. Comparison 
with the gia nt and dwarfs spectr a from the IRTF spec¬ 
tral library (Rayner et al. 2009), reveal that they are 
most probably K-type giaiit stars, although in the case 
of #33 (16114657-3906047) and #34 (16120495-3906256) 
we cannot with certainty exclude that they are young, 
since the CO bands are not as deep as in the other two 
objects. Additionally, #33 exhibits Bry in emission, i.e. 
this object might be a stellar member of Lupus 3. For 
the remaining objects listed in Table [^without any com¬ 
ment, the low resolution and S/N of our spectra prevents 
a straightforward classification. 


3.2. Model fitting 

To determine the effective temperature and extinction 
of the three objects showing the prominent H-ban d peak, 
we fitted t he spectra with the BT-Settl models (Allard 
et al.|2011|). The procedure is identical to the one applied 
in other papers of the SONYC series using NIR spec- 
troscopy (for a deta iled description of the procedure, see 
Scholz et al.|2012b). The best fit solution is searched on 
a grid of 'Jeff varied between 2000 and 4000 K in steps of 
100 K, and Ay within ±5 mag of the value derived from 
photometry, by assuming the intrinsic (J — K)q = 1, and 


extinction law from Cardelli et al. (1989 1“ The logg 


is kept at a single value of 3.5, suitable tor young ob¬ 
jects that are still contracting. For the two objects in 
Lupus 3 no, or a very small adjustment to the photo¬ 
metric value of extinction (AAy < 1) is needed to pro¬ 
duce a satisfying fit. For SONYC-ChaI-1, on the other 
hand, the value about 2 — 2.5 mag lower provides a much 
better solution. For the uncertainties of the fitting pro- 


As argued in|Scholz et al. (|2009j, the value (J — K)o = 1 is 
appropriate for objects ot the M spectral typ e. This is in agree¬ 
ment w ith values for 5-30 Myr old M-dwarfs by|Pecaut &; Mamajekl 
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ISO-Chal-192 


H,0 


CO 


SONYC-Lup3-29 2800K Av=3.5 




SONYC-Lup3-28 2800KAv=1.0 

■ - 



ol_ I _^_L 


Fig. 2.— Spectra of the confirmed members of Cha-I and Lu¬ 
pus 3 (black), along with the BT-Settl best-fit models shown in 
red. All the spectra are normalized at 1.6^m, and shifted by ar¬ 
bitrary offsets for clarity. The object spectra are as observed, and 
the models have been artificially reddened. The region strongly af¬ 
fected by telluric absorption (1.8 — 1.95 fim) is not shown for clarity. 
The prominent molecular absorption bands observed in ultracool 
objects are marked. The lines marked at the top of the plot are 
the H 2 ro-vibrational transitions. 


10"* 


10''® 


< 10*" 
"e 

10''* 


10"'* 


10'" 

1 10 
X [^m] 

Fig. 3.— Spectral energy distribution for two objects from our 
follow-up showing infrared excess. The solid line in the left panel 
is a BT-Settl atmospheric model for a low-gravity object with 

Teff=2800K. 


ISO-ChQl-192 


oSONYC 
•2MASS 
OWISE 
■ Spitze'' 



cedure, we adopt the conservative values of ±200 K for 
the Teff, and ±lmag for the Ay, same as in our ear¬ 
lier works. SONYC-Lup3-28 was previously confirmed 
as a p r obabk substellar member of Lupus 3 by |Comer6n| 
et al.j (|2013D, who derived the = 2800K, in agree¬ 
ment wTEH~the value derived here. SONYC-ChaI-1 was 


previously reported by Luhman (20071, and classified as 
a probable member of Cha-1, with spectral type >M9, 
and Teff <2400 K, in agreement with our results. 

3.3. Spectral types 

We derived spectral types for the three late type mem¬ 
bers confirmed in this work, usi ng the following spe ctral 
indices: (1) H-peak index (HPI , Scholz et al.|2012b ), (2) 
H 2 O index (lAl lers et al.||2007|), and (3) (j-index (iVVilk- 
ing et al.|l9y9). 'i'he three indices give consistent JNiK 
spectral types for the two M-dwarfs found in Lupus 3. 
For SONYC-Lup3-28 we get M8 (HPI and Q-in dex), and 
M7 (H 2 O), i n agreement with M7.5 derived by|Comer6n 
et al.| ( 2013|) . For SONYC-Lup3-29 we get M8 (HPI), 
Alb (H 2 U), Ld M7 (Q). For SONYC-ChaI-1, the indices 
result in L2 (H 2 O), and L4 (HPI and Q-index). While 
all the t hree indices are defi ned for young BDs earlier 
than LO, Scholz et al. (2012 b|) note that HPI mi ght hold 
also at later spectral types. jAllers et al.j (|2007|) derived 
two SPT-H 2 O index relations, one tor young BDs in the 
AI5 - LO range, and other for the field dwarfs in M5 - 
L5 SpT range. They note that the two relations appear 
very similar, and thus the SpT-index rel ation can also be 
used for the young BDs later than LO. Luhman (20071 
derive spectral type >AI9 for SONYC-Chal-1. 

Prior to the calculation of the spectral types, the spec¬ 
tra had to be corrected for extinction. To determine the 
influence of extinction on the resulting spectral types, 
we varied the Ay derived in the previous section by 
±lmag. For HPI, we obtain the uncertainty of ±0.3 
subtypes for SONYC-Lup3-28 and -29, and ±0.5 for the 
noisier SONYC-ChaI-1. For H 2 O index, we get ±0.3 for 
all three objects, while the Q-index does not depend on 
the extinction. Additionally, the uncertainties for the 
spectral types derived from the three indices from their 
respective defining papers are ±0.4 for HPI, ±1.0 for the 
H 2 0-index, and ±1.5 for the Q-index. These uncertain¬ 
ties have been added in quadrature to the uncertainties 
from the extinction, and taken into account for the calcu¬ 
lation of the final spectral types, which were rounded-up 
to the nearest half-spectral type. 

In Figure we show a comparison of the SONYC- 
ChaI-1 spec trum with the spec tra of field ultracool 
dwarfs from jAllers & Liu (2013), and wi th the young 
L-type members in Upper Sco and pOph. Allers & Liu 
(2013) presented a spectroscopic study of field ultracool 
dwarfs having spectroscopic and/or kinematic evidence 
of youth (« 10 — 300 Myr). The objects are divided in 
three gravity classes based on the shape and strength of 
various features in their spectra. Three L3 objects ex¬ 
hibiting high-, intermediate-, and very-low gravity fea¬ 
tures are shown in the left panel of Figure At the top 
of the same panel we also show a very-low gravity L4 ob¬ 
ject to facilitate the comparison with the young members 
of star forming regions, shown in the right panel. There 
are no s uitable intermediate -gravity L4 standards pre¬ 
sented inlAllers & Liu (2013). From the shape of the H- 
eak" 


band peak, it is evident that 


SONYC-Chal-l cannot be a 
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16100597-3908062 




16103359-3908214 


'lOO-iA ' 


16101171-3909198 


16094266-3903207 



16095268-3906404 ' 


16094374-3907119 


16094720-3908424 


T~~~^ 


16082587-391 


16080396-391 



1.8 2.0 2.2 2.4 



1.6 1.8 2.0 2.2 2.4 


Fig. 4.— Observed spectra of candidate objects in the direction of Lupus 3, for which the membership cannot be determined with 
certainty. While it is clear that none of these objects is a young substellar member, it cannot be excluded that some of these still could be 
stellar members with spectral types earlier than ~M5. All the spectra are normalized at 1.6/im, and shifted by arbitrary offsets for clarity. 
Three spectra shown at the bottom of the left panel have lower S/N and were smoothed to better appreciate the overall broad shape of the 
spectra. The region strongly affected by telluric absorption (1.8 - 1.95 //m) is not shown for clarity. 


11062804-7735547 





I : 1 

o 3- I 


2r 


11083077-7745505 


11085038-7745530 


1.6 1.8 2.0 2.2 2.4 



Fig. 5. — Observed spectra of candidate objects in the direction of Cha-I, for which the membership cannot be determined with certainty. 
While it is clear that none of these objects is a young substellar member, it cannot be excluded that some of these still could be stellar 
members with spectral types earlier than '*^M5. All the spectra are normalized at 1.6/i.m, and shifted by arbitrary offsets for clarity. With 
four exceptions, the spectra were smoothed to a lower resolution to suppress the noise. The region strongly affected by telluric absorption 
(1.8 - 1.95 /im) is not shown for clarity. 
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with the spectra of field ultracool dwarfs from|Allers Liu]<|2013 
left panel), and young L-type members of star forming regions 
(right panel). Spectra of the comparison objects are shown with 
no reddening (grey), and with the extinction of Ay = 2.5 (the 
value estimated for SONYC-ChaI-1 from the model fitting; ma¬ 
genta) applied to it. The field objects are, from bottom to top, 
2MASSJ0407-hl546 (L3 field gravity), 2MASSJ1726-hl538 (L3 in¬ 
termediate gravity), 2MASSJ1615-I-4953 (L3 very-low gravity), and 
2MASSJ15515237-h0941148 (L4 very-low gravity). The young ob- 
jects include two me mbers of Upper Sco, UScoJ16 0603-221930 (L2: 
Lodie u et al.| 2008ll and 1RXSJ1609-2105B fL4: iLafreni ere et al. 
2UiU|)' and U W IR-Oph33, a member of pOph fL4; jAlves de 
Uliveira et al. ] |2012| | . For a reference, we also show the best-ht 
hi i-bettl low-gravity model (2200 K). All spectra are smoothed to 
the same resolution, and normalized at 1.7pm. 

normal field dwarf, but it is difficult to judge whether an 
intermediate or a low-gravity atmosphere provides a bet¬ 
ter fit. Overall, the spectral features match well both L3 
or L4 low-gravity objects, but the slopes between H and 
K bands are different, in the sense that our object ap¬ 
pears bluer than any of the field templates, even without 
extinction. 

To date, there is only a limited sample of very young 
L-type members of star forming regions with good qual¬ 
ity spectra available for comparison, especially at spec¬ 
tral types later than LI. In the right panel of Figure 
we show two members of the Upper Scorpius ass ocia- 

and 


Lodieu et al. 


tion, UScoJ160603-221930 (L2; . 

1RXSJ1609-2105B (L4; ILafreniere et al.||2LilUD, along 
with CFHTWIR-Oph33, an L4 member of p Upn (Alves 


de Oliveira et al. |2012 ). The H-band portion of the 
SUJN Y (J-Chal-1 spectrum matches well the young L4 
templates, but in the K-band it again appears to be more 
consistent with Ay = 0, rather than with Ay = 2.5, the 
value estimated from the spectral model fitting. For a 
comparison, on top of the right panel we show the best- 
fit BT-Settl model, which clearly prefers the higher value 
of the extinction. 

We note that detailed spectral type comparisons at 
young ages can be complicated by possible excess emis¬ 
sion from disks o r accretion. As demonstrated by |Daw- 


son et al. (2014), these characteristics alter the spectra 
of Class ii objects, which then show different normalized 
levels between individual NIR bands. The effect is also 
present in Class III object spectra, to a somewhat lesser 


extent. 

For now, we conclude that the spectrum of SONYC- 
ChaI-1 is consistent with a spectral type L3-L4, Ay < 3, 
and a low-gravity atmosphere. More high quality spectra 
of young L-type objects are needed in order to establish 
a proper NIR spectral type sequence, and to create a 
gravity sequence at each spectral type, in c omparison to 
the field ultracool low-gravity dwarfs from|Allers & Liu 

(plk. 


3.4. Hertzsprung-Russ ell Diagram 

In Figure we show HR diagrams for Cha-I (left 
panel), and Lupus 3 (right panel). The luminosities were 
calculated from the J-band magnitudes corrected for ex¬ 
tinction, distance modulus, and bolometric correction. 
We adopt a distance of 200 pc for Lupus 3, and 160 pc 
for Cha-I. Extinction was calculated from the J — K col¬ 
ors, assuming the intrinsic colors according to spectral 
type o f each member, as listed in Pecaut fc Mamajek| 
(2013) for young ( 5-30 Myr) dwarf stars. We use the 


extinction law from Cardelli et al. 
The bolometric correction in the ■ 


(198^, with Ry = 4. 

3 


band (BCj), for the 
objects with > 2750 K, is calculated from t he poly¬ 
nomial relation bet ween BCj and Teg derived by |Pecaut] 
& Mamajekj (|2013|). For the objects with Teg < 2750 K, 
where this relation is not valid, we determine the BCj 
as a function of J — Ks color, as derived in Schmidt 


et al. (2014). The former relation is valid for young (5- 
30 Myr) dwarfs with spectral type earlier than M6, while 
the latter was derived for field dwarfs with spectral types 
M7-L8. 

The HR diagram for Lupus 3 i s an u pdated version of 
the one presented in|Muzic et al.|(|2014|), where we added 
the two objects confirmed in this work (red diamonds), 
and the candidate member J16114657-3906047 (purple 
square). For the latter we can only infer a lower limit on 
the Teg. Since we do not know the exact spectral type 
of this object, we adopt the intrinsic J — K color 0.6- 
0.9 (used to calculate the Ay), and J-band bolometric 
correc tion 1.5 ± 0.2 mag, both s uitable for young K-type 
stars ( Pecaut fc Mamajek|2013 ). To construct the Cha-I 
HR diagram, we us e the photomet ry, Tgg, and extinction 
from the census by Luh man j 2007), add the o bjects from 
Luhman et al. (2008~ Luhman & Muench| (2008), and 
the two sources contirraed in this work (red diamonds). 
The typical error-bars are shown in the lower left corner 
of the two plots. Teg estimates for Lupus 3 come from 
various works, and we take ±200K as a representative 
uncertainty (see Muzic et al. 2014 for details). The Tgg 
of the Cha-I members marked with crosses comes from 
the Luhman census, who converts spectra l types to Tgg 
using the scale from Kenyon & Hartmann (1995) for the 
stars with spectral ^pes e arlier than M, and the one 
from Luhman et al. (2003) for the M-type. The latter 
scale was des igned to be compat ible with the evolution¬ 
ary models of Baraffe et al. (1998). Therefore, in addition 
to the uncertainties m spectral types, these temperature 
estimates are subject to a systematic uncertainty in the 
temperature scale, which is difficult to predict, since it 
depends on vari ous details of the models, a nd is probably 
at least ±100 K (Luhman & Muench|2008 ). We therefore 
decide to adopt a single value of ±200 K tor all objects in 
the plot. Different uncertainties in luminosity for objects 
hotter, or cooler than 2750 K in Figure]^ result from the 
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TABLE 3 

Confirmed members of Cha-I and Lupus 3, from the SofI follow-up. 


No.‘‘ 

name 

a(J2000) 

5(J2000) 

SpTb 

Tefl/K 

Au/mag 

mass/MQ^^ 

log(L/L 0 ) 

comments 

6 

SONYC-Chal-l 

11 08 30.31 

-77 31 38.6 

L3 ± 0.5 

2200 

2.5 

0.009-0.012 

-3.44 ±0.10 

>M9. <2400K'i 

11 

ISO-ChaI-192 

11 09 28.57 

-76 33 28.1 





- 0.11 ± 0.20 


19 

SONYC-Lup3-28 

16 08 33.04 

-38 52 22.7 

M8±0.5 

2800 

1.0 

0.05-0.06 

-l.37T0.il 

M7.5, 2800 K® 

24 

SONYC-Lup3-29 

16 10 01.33 

-39 06 45.1 

M7.5 ± 0.5 

2800 

3.5 

0.05-0.06 

-0.81 ± 0.11 

IR excess 


^ Same as in Table 
NIR spectral type 


estimate based 

on Teft 

|Luhman| J2007 


IComeron et al. 

J2013ll 


uncertainties in BC.i that were adopted from two dif- 
ferent works (Pecaut & Mamajek 2013| Schmidt et al. 
2014| . The two sources with uncertain spectral type 


have slightly larger logL uncertainty, as a consequence 
of a range of extinctions and bolometric magnitudes. 

In both plots, there are sources appearing well below 
the main sequence formed by the cluster members. In 
case of Cha-I, for 7 out of the 8 sources with Teg > 3000 K 
that are loc ated below the 100 Myr isochrone, |Luhman] 


(2004 


20071 argues in favor of these sources being gen- 
uine members of Cha-I, based on mid-IR excess and/or 
strong emission lines. Only one of the 8 sources, OTS 

32, has an uncertain membership. In Lupus 3, 3 of the __ 

underluminous objects also show presence of a disk and Tpff < 2800K (Helling et al. 2008 Helling fc Casewdl[ 


man et al.| (|2011a|b|) , studying two planetary-mass com- 
panions to young stars, suggest that the effect might be 
caused by an overestimate of the They were able to 
reproduce the observed properties of the two objects by 
modeling low-gravity, cloudy atmospheres experiencing 
non-equilibrium chemistry, and requiring the Teg that is 
several hundred K lower than reported earlier. Therefore, 
the position of the coolest objects in the HR-diagram 
with respect to the theoretical isochrones can reflect a 
problem with the bolometric luminosities and/or the at¬ 
mosphere models in this Teff regime. The latter is possi¬ 
bly due to treatment of dust in various models, as clouds 
start to form i n the atmospheres of cool dwarfs with 


H„ emission, but no mid-IR excess (Muzic et al. 

2014 

Comeron et al. 

2003 Fernandez & Comeron 

2005 

iTM 

discussed in 

Comeron et al. (2(J(J3|), the underluminosity 


20141. 


plained by the object being seen only through scattered 
light (edge-on disk, or embedded Class I sources), or by 
the accretion-modified evolution, similar to what is de¬ 


scribed in Hartmann (19981 and Baraffe et al. (2009). 
The examples of the detail modeling ot young uiiderlu- 
minous objects with edge-on disks can b e found in e.g. 


Scholz et al. (2010); Huelamo et al. (2010); Petr-Gotzens 

let al.|(|20i0| - - - - 

In Cha-I, the sequence of the objects hotter than 
^ 2700 K is following the isochrones, and is mostly con¬ 
fined between those marking 1 and 10 Myr. Most of 
the coolest objects, on the other hand, fall between 10 
and 100 Myr, i.e. they appear systematically underlu¬ 
minous relative to theoretical isochrones. One possible 
culprit are the bolometric corrections used to calculate 
the luminosities. To cal culate BC.i for these objects, 
we use the relation from Schmidt et al. (2014), which 
was derive d for the field d warfs. As summarized pre¬ 
viously in |Luhman (2012), the bolometric corrections 
for young L and i' dwarfs are probably not the same 
as those for standard field dwarfs, which may lead to 
underestimated luminosities for the young dwarfs. The 
bolometric corrections appropriate for the young late M 
and L dwarfs are not yet available, but a few indi vid- 
ual values do exist in the literature. For example, |Za-| 

or 


patero Osorio et al. (2010) derive BCj=1.16 ± 0.10 
a young L3 dwarf G196-3B, whereas a field dw arf of the 


same spectr al type would have BCj=1.85±0.14 (Schmidt, 
et al.| 2014). This difference in BCj shifts an object by 
~ 0.3dex in the positive y-direction, i.e. enough t o fall 
above the 10 Myr isochrone. On the other hand, |Bar-| 


4. DISCUSSION 
4.1. Chamaeleon-I 

Based on the survey presented in this work, we can 
put limits on the number of VLM sources that are 
missing in the current census of YSOs in Cha-I. As in 
other works from the SONYC series, we use the suc¬ 
cess rates of our spectroscopic follow up (defined as a 
number of confirmed members divided by the number of 
the spectroscopically surveyed candidates), to estimate 
the expected number of members among the remain¬ 
ing photometric candidates. The listed uncertainties are 
the 95% confidence intervals (Cl), calculated using the 
Clopper-Pearson method, which is suitable for small- 
number events, a nd returns conservative CIs comp ared 
to other methods ( |Gehrels||l986 Brown et al.||2001 |. 

While until now for this estimate we always considered 
the entire candidate lists above the completeness limit, 
in case of Cha-I it might be more suitable to divide the 
candidate list in two, and consider the candidates with 
I < 21, and those with 21 < J < 23 separately. The 
reason for this is the following. The candidate list with 
/ < 21 contains 46 objects. We took 16 spectra with 
VIMOS, and confirmed 7 members. This gives a success 
rate of 7/16 = 44+|4% (95% CIs). Among the candi¬ 
dates with 21 < J < 23 that were specifically targeted 
in this most recent follow-up, we took in total 17 spectra 
(15 with SofI and 2 with VIMOS), and confirmed only 
I probable VLM members of Cha-I (the other red dia¬ 
mond in Figure [l] marks ISO-ChaI-192, whose spectral 
type is unclear). The success rate is therefore 6tlg^%, 
which is significantly lower than for the candidates with 
/ < 21. The difference cannot be caused simply by dif¬ 
ferent wavelength domains used in spectroscopy (optical 
versus NIR), as all but one of the members confirmed 
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lOMyr.... 
.. ,100Myr 


5Mxr 
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SONYC confirmed ♦ 
J16114657-3906047 ■ 


. X ^ 


IS X jfx' 


: I- 


H Teff< 2750 K 


-3 r I- 1 - 1 Teff >2750 


T.JK] 


4000 3500 

T.„[K] 


Fig. 7. — H-R diaeram for Cha-I (left) and Lupus 3 (right). The crosse s mark the member s of Chal from|Luhman|l|2007|l;|Luhman et al.| 
l|2008|l; [Luhman &: Muench| ||2008[l and Lupus 3 from the census table in [Muzic et al.| (|2014|l. Red diamonds mark the members conhrmed 


in tms work, and tne purple square marks the candidate member of Lupus 3 ibii4b5Y-3yUbU47. The arrows signal the lower or upper limits 
on Tggf. The dashed lines show BT-Settl theoretical tracks for ages of 1, 5, 10, and 100 Myrs. 


with VIMOS have Teff below 3500 K, and would have 
also been recognized in the NIR. 

From the success rates mentioned above, we can now 
estimate the number of missing VLM members in our 
candidate list. The number of unconfirmed members 
among the remaining candidates in a sample can be cal¬ 
culated by subtracting the number of spectra taken by 
SONYC and number of objects observed by other groups 
from the total number of candidates in the sample, mul¬ 
tiplied by the success rate. In the higher brightness bin 
(/ < 21), among our 46 candidates without SONYC 
spectra, we find additional 8 obj ects class i fied a s mem - 


bers, and 4 non-members from Luhman (2004 2007). 


Thus, the number of missing VLM objects with [I <21) 
is (46 — 16 — 8 — 4) X 0.44jlg;24 = 8 + 4 . If we now consider 
the 21 < / < 23 part of the Iz CMD, it contains 60 can¬ 
didates, of which 17 have SONYC spectra. One probable 
VLM member was confirmed by SONYC, and one addi¬ 
tional VLM was confirmed in other surveys. We estimate 
therefore the expected number of members in this bin to 
be (60 — 17 — 1) X 0.06l(Q;g| = 3+3°. We expect most of 
these missing members to be substellar, although some 
of them might be more massive, embedded members of 
Cha-I. 

With the above estimates in mind, we can say that the 
census of VLMOs in Cha-I is almost complete, down to 
0.004 — 0.009 M©, for Ay <5 and age of 2 Myrs. 

4.2. Lupus 3 

Following the dis cussion in the Sections 5.3 to 5.5 of 


Muzic et al. (|2014), we can update the mass function 
based on the lollow-up presented in this work. For this 
we take into account only the candidates above the com¬ 
pleteness levels of our survey. The updated success rate 
in the “/J-pm” sample is therefore 9/50 = 18“'l9^%, while 
for the photometric-only candidates we have 0/81 = 
Gtl%. As in the previous section, the number of un¬ 
confirmed members among the remaining candidates in 
a sample can be calculated by subtracting the number 
of spectra taken by SONYC and number of objects con¬ 
firmed by other groups from the total number of can¬ 
didates in the sample, multiplied by the success rate. 
Therefore, the estimated number of unconfirmed mem¬ 
bers in our survey is (53 - 50) x 0.18lj];J| -F (337 - 53 - 
81 — 1) X 0“*'°'g^ « 1~L\°. The first term in this equa¬ 
tion refers to the “/J-pm” sample, the second to other 


(photometric only) candidates, and the quoted error is 
the 95% CL We can repeat the same calculation to get 
an approximate number of missing BDs above the com¬ 
pleteness limit, but the result is essentially the same as 
for all VLM objects, i.e. any of the 0-11 missing objects 
could be either stellar or substellar. 

Based on these numbers, it is evident that very few 
potential VLM members of Lupus 3 are left to be dis¬ 
covered, at least down to the completeness limit of our 
survey, which is equivalent to 0.009—0.02 Mq, for Ay <5. 
Our survey covers a large portion of the entire Lupus 3 
cloud, and encompasses the entire high-extinction band 
around the two H erbig Ae/Be members (see Figure 1 in 
Muzic et al.|2014|) where majority of the known members 
are located. 

Comparing Lupus 3 with Cha-I, we note that the suc¬ 
cess rate of spectroscopic confirmation of VLM candi¬ 
dates is significantly lower in the former, despite hav¬ 
ing used a homogenous method for selecting candidates. 
This can be explained by different contamination rates 
by background objects along the two lines of sight. The 
galactic longitude of Lupus 3 is Z ~ 340°, which means 
that we are looking towards the inner, more densely pop¬ 
ulated parts of the Galaxy, compared to regions along 
the line of sight to Chamaeleon (Z ^ 297°). The Be- 
sancon Milky Way stellar population synthesis model 
(Robin et al. 2003) yields about five times more objects 
in the direction of Lupus 3, than towards Cha-I, within 
the same area on the sky and for the same photometry 
limits. Lupus 3 is at the same time less populated: the 
current stellar census in this region contains roughly 80 
members, while in Cha-I there are ^180. Higher contam¬ 
ination by background sources, combined with a smaller 
overall population is therefore the most plausible cause 
for different success rates in our spectroscopic surveys. 

4.3. Distribution of spectral types 

In Figure we show the distributions of spectral types 
for the spectroscopically confirmed VLM population in 
Lupus3 (left panel), and Cha-I (right panel), found 
within the region covered by our surveys. The solid er¬ 
ror bars are Poissonian confidence intervals calculated 


with the method described in Gehrels (1986). We do not 
try to correct the histograms tor the missing objects es¬ 
timated in the previous section, because it is not clear 
which spectral type they would have. Assuming that 












































SONYC IX: Cha-I & Lupus 3 


11 



25 


20 - 


15 - 


10 - 


5 - 


Cha-I 0.5M®T 0.08M® 0.02M® 0.006M® 


I 


I 


I 


I I I I 


I 




MO M2 M4 M6 M8 LO L2 L4 
Spectral type 


Fig. 8.— Distribution of spectral types for the VLM population of Lupus3 (left) and Cha-I (right), within the area of our surveys. The 
rough mass limits according to BT-Settl isochrones are shown on top of the plot; the dashed lines mark the completeness limit of the two 
surveys for ArA=5. 


they have moderate extinction (Ay < 5), we could dis¬ 
tribute them in the substellar bins of the two histograms. 
However, our previous surveys have revealed a significant 
number of more embedded stellar members in the same 
part of the CMD where those moderately-extincted sub- 
stellar members are found. In any case, the numbers of 
missing members are low, especially in Lupus 3, and they 
cannot strongly affect the spectral type distribution. 

The overall trend seen in the two histograms is very 
similar, with a rise in the number of objects at spectral 
types earlier than M5, a peak at M5, and a drop in the 
number of objects at later spectral types. A similar be¬ 
haviour is_obsCTyed_hi other st ar forming regions, e.g. 
IC 348 (|Luhman||2007|) , p Oph (|Alves de Oliveira et al. 
' 2012| , and NGC 1333 ( [Scholz et al.|2012a| |"^" The sharp 
drop in the number of objects at spectral types M7 and 
later is seen in all these distributions, and is certainly a 
real feature of the IMF, and cannot be attributed to the 
incompleteness of each survey, especially in Cha-I where 
we are complete down to ^L3. The survey in Lupus 3 is 
somewhat shallower, but from the analysis in Section [4^ 
it is evident that there are very few new members left 
to be discovered in this region down to ~L0. Even if 
we assume that all the objects missed by our surveys 
are substellar (i.e. in the bins M7 and later), the drop 
at late spectral types would remain. The two deepest 
SONYC surveys, in NCC 1333 and Cha-I, both contain 
very few objects later than M9. The two surveys are 
both complete down to ~0.005Mq, and reveal that the 
free-floating objects below D-burning li mit a re rare. This 
will be further discussed in the Section mH 


4.4. Star/BD ratio 

To assess the numbers of stars and BDs, we use the 


approach described in Scholz et al. ( 2013 ). In short, by 
comparing the multi-band photometrj^ with the pre¬ 
dictions of the BT-Settl evolutionary models, we derive 
best-fit mass and Ay for each object, for the assumed dis¬ 
tance of 160 pc and age of 2 Myr for Cha-I, and 200 pc 
and 1 Myr for Lupus 3, and the extinction law from 


See Figure 13 in 


Muzic et al. 


1 2014| for a direct comparison. 


IJHK photometry complemented with the optical BVR 
bands where available from public catalogs. 


TABLE 4 

Star/BD number ratio 


Chamaeleon-I 

stars 

BDs 

RC 

R2‘= 

all 

< 1M0 

< 1M0 

all 

all 

> O.O3M0 

2.4 - 3.8 
2.0 - 3.4 
3.0 - 6.1 

1.7 - 2.8 

1.5 - 2.5 

2.5 - 5.6 

Lupus 3 

stars 

BDs 

RA 

R2‘= 

all 

< 1M0 

< 1M0 

all 

all 

> O.O3M0 

2.3 - 3.4 
1.7- 2.7 
3.0 - 6.1 

2.0 - 3.0 
1.6 - 2.5 
2.9 - 5.8 


^ in the area covered by the Luhman census 

^ in the area of our survey _ 

^ in th e ar ea covered by surveys o f |Merin et al.| 
j2008^ and [Comeron et ah] l|2009^ 


Cardelli et al. (19891 with Ry = 4. 

For the stellar-substellar mass boundary, we take the 
value at the solar metallicity, O.O75M0. All the objects 
with estimated masses below 0.065 Mq are counted as 
BDs, and with masses above 0.085 Mq as stars. The 
remaining objects at the border of the substellar regime 
are once included in the higher mass bin, and then in the 
lower mass bin, thus resulting in a lower and upper limit 
of the star/BD ratio. 


20071), and complement it with the obiects identified in 

Luhman et al. 

(2008): Luhman & Muench 

(12008), re- 

solved binaries trom 

Daemgen et al. (21)131, and one VLM 


number ratios, calculated for different mass bins, and for 
the area of the Luhman census, as well as the area en¬ 
compassed by our survey. The limits on stellar mass 
< IMq, and BD mass > O.O3M0 were chosen t o facil¬ 


itate comp arison with our previous works (e.g. Scholz 


i'he stellar and substellar population of Lupus 3 has 
not been as thoroughly surveyed by spectroscopy in the 
past as that of Cha-I. If we look at only the spectroscop- 
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ically confirmed members, it is clear from Fig. [^that the 
number of objects in the M6 bin is larger than the to¬ 
tal number of objects in lower mass bins. Since a good 
fraction of the objects in the M6 bin are at the substellar 
border, we end up with a very large span of values for 
the star-to-BD ratio. However, Lupus 3 certainly has a 
larger fraction of candidate sources that still await a spec¬ 
troscopic con firmation. For example, the work of |Merm 


et al. (2008) contains 124 candidates, of which only 46 


have been surveyed by spectroscopy. We therefore ap¬ 
ply a slightly different approach to calculate the star-to- 
BD ratio in Lupus 3. We use the following (candidate) 
member lists: (1) ce nsus of spectroscopi cally confirmed 
M-type objects from Muzic et al. (2014), updated with 
the newly discovered object SU M Y (J-Lup3-29; (2) candi¬ 
dates from Men'n et al. (2008), excluding those already 
in the census table. I'he numbers from this sample are 
multiplied by the success rate o f the spectroscopic follow¬ 
up (36/46; Mortier et al.|20ir ), and by the factor 1.25 to 
account for the tact that the survey is only sensitive to 
objects with disks, whose fraction is estirnated to be 70- 
80%; (3) candidates from Comeron et al. (2009), exclud¬ 
ing those already in the census table. I'he numbers from 
this sample are multiplied by 0.5, which is the success 
rate of the s pectroscopic follow-up presented in|Comer6n| 


et al. (2013); (4) kno wn membe r s with the spectral type 
earlier than M, from Comeron (2008). The results for 
the stars-to-BD ratio are shown in 'Fable for the area 
encompass ed by our survey, a s we ll as for the larger are a 
covered by|Merin et al.| (|2008|) andjComeron et al.|(|2009|) . 

The nuinbers listed m Table tor Cha-1 and Lupus 3 
are consistent with each other, and are also generally con¬ 
sistent with th e star-to-BD ratios d erived for NGC 1333 
and 1C 348 by Scholz et al. (2013). Clearly, there are a 
number of uncertainties involved in this calculation, such 
as the choice of the isochrones and extinction law used to 
derive masses, or uncertainties in distances to star form¬ 
ing regions. In Lupus 3, in particular, there could also 
be some effects of possible incompleteness at the overlap 
of the different studies. 

The numbers in Table |4] do not include our estimate for 
the objects that are “missing” in our census. In Cha-I, 
where the stellar population has been thoroughly char¬ 
acterized by spectroscopy in the works of Luhman et al., 
we estimate that there might still be IFts^ unidentified 
members in our selection box. According to the BT-Settl 
models, at moderate extinctions (Ay < 5), these objects 
should be substellar, although some of them might be 
more embedded stellar members. If these objects would 
be substellar, the ratio of the star-to-BD ratio R 2 could 
be as low as 1, with the upper limit of 2.3. In Lupus 3, 
where we surveyed almost entire population of the pho¬ 
tometric and proper motion candidates, the number of 
missing objects is With this upper 95% confidence 
limit, the star-to-BD ratio R 2 would be slightly lower: 
1.7-2.4 taking into account all stars and BDs, and 1.4- 
1.9 if we consider only stars below IM©. 

To conclude, the star/BD ratio in Cha-I and Lupus 3 
is consistent with the findings in other clusters, and we 
find that for each formed BD there are between 1.5 and 
6 formed stars. To probe finer differences, one needs 
a more sophisticated analysis and better constraints on 
fundamental parameters. 



Fig. 9.— Mass spectrum for low mass stars and BDs in Cha-I 
(upper panel) and Lupus 3 (lower panel), within the areas of our 
surveys. The bin sizes were chosen to contain similar number of 
members. The red squares in the Cha-I plot account for the num¬ 
bers of missing objects (see text). The vertical error bars represent 
the Poisson uncertainties. The solid black lines show the best fit to 
the data (black circles), and the red dashed line in the upper panel 
is the best-fit to the data when the red squares are used instead 
of the circles at their respective masses. The vertical dotted lines 
mark the approximate locations of H- and D-burning limits, and 
the dash-dotted lines represent the completeness of the two surveys 
at Ay = 5. The grey dashed line in the lower panel is the Kroupa 
segmented power-law mass function. 



mass (M©) 

Fig. 10.— Mass spectrum in Lupus 3. Black circles and his¬ 
togram show the same data as shown in the lower panel of Figure]^ 
only restricted to masses above O.O2M0, which is our conservative 
completeness limit at Ay = 5. Red circles show the mass spectrum 
for the same sample of stars, but with masses calculated after ran¬ 
dom variation of Av^ (see text). Blue crosses and histogram show 
the mass spectrum after adding 11 substellar objects with masses 
above our completeness limit at Av^=0 (0.009 Mq). The bins are 
modified to have similar error bars, which are not shown here for 
clarity, but are similar to those in Figure]^ 


4.5. Initial Mass Function 

With the masses calculated as described in the pre¬ 
vious section, we can construct the IMFs in Cha-I and 
Lupus 3. They are shown in Figure where we varied 
the binsize in order to achieve similar numbers of mem¬ 
bers in each bin. We choose to express the IMF in the 
power-law form, as dN/dM oc M“". The data needed 
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to reproduce the plot are given in Table In Cha-I the 
best fit yields a = 0.68 ± 0.09 (black solid line). We 
can also account for the number of obj ects missing in 
our census, as estimated in Section |4.l| It is not sim¬ 
ple to predict the masses that these “missing” objects 
should have, because of the variable extinction in star 
forming regions, where more massive embedded stellar 
members might mix in the color-magnitude space with 
the less embedded substellar members. However, if we 
assume the typical Ay < 5, then our missing members 
should be substellar. The objects in the 21 < / < 23 bin 
were included in the lowest-mass bin (0.001 — O.O 3 M 0 ), 
as they are expected to have masses below 0.015 Mq. 
Objects with 17.5 < / < 21 should have spectral type 
M7 - LO, and we distribute half of those in each of the 
two lowest-mass (substellar) bins. With this correction, 
the best-fit power-law index (dashed red line) becomes 
a = 0.78±0.08. The slope remains essentially unchanged 
if we put all the missing objects with 17.5 < / < 21 to the 
more massive of the two substellar bins (O.O3-O.O75M0). 
Despite the caveat mentioned above, this short exercise 
tells us that even in the (somewhat extreme) case that 
all the missing members were substellar, the slope of the 
IMF would not be strongly affected. 

In Lupus 3, we get the best fit single power-law in¬ 
dex a = 0.79 ± 0.13 (bottom panel of Figure]^. Here 
the lowest-mass bin contains all the substellar objects. 
While in Cha-1 all the data points seem to be consistent 
with a single power-law, in Lupus 3 there might be an 
indication of the flattening of the slope in the substel¬ 
lar regime. For a compari s on, w e show the segmented 
power-law IMF by Kroupa (2001), where the power-law 
slope changes at 0.5 and 0.08 M 0 . The Kroupa IMF, 
normalized to match the middle point of the Lupus 3 
IMF at O.I 55 M 0 , is shown by the dashed grey line in 
the lower panel of Figure with the slopes a = 2.3 
for M> O. 5 M 0 , a = 1.3 for 0.08 < M < O. 5 M 0 , and 
a = 0.3 for M< O.O 8 M 0 . The possible turnover in the 
power-law mass function in the substellar regime that 
we observe in Lupus 3, has not been, to our best knowl¬ 
edge, reported in any other spectroscopically confirmed 
sample of a young cluster. One possible exception is the 
young, sparsely populated association ry Cha, which be¬ 
side the population of ~15 stars of spectral ty pe K and 
M, does not con tain any BDs down to O.O15M0(Luhman 
fc Steeghs|2004 ). Assuming the star-to-BD ratio of 1.5 
6 that we hnd for Lupus 3 and Cha-I, one would expect 
a population of 3 - 10 BDs. However, a comparison of 
the mass distributions in 77 Cha, Cha-I, and IC 348 by 
means of a two-sided Kolmogorov-Smirnov test, does not 
reveal a significant d ifference between the three clusters 
( Luhman et al.|[2009 ). 


to contain similar numbers of objects. Black circles and 
histogram show the spectroscopic IMF from Figure]^ re¬ 
stricted to masses above our conservative completeness 
limit at Ay = 5 (O.O 2 M 0 ). We see that even after adding 
the upper limit of the possibly missing objects, the flat¬ 
tening of the slope remains. 

A second suitable test concerns possible influence of the 
extinction. Here, we derived masses and Ay simultane¬ 
ously by fitting the photometry to the predictions of the 
evolutionary models. Alternatively, we could have fixed 
the Ay to the value derived from spectroscop y, or from a 
single intrinsic color (e.g. J— AT, as in Section 3.4). Since 


for the majority of the objects these values do not differ 
by a large amount, the resulting IMF looks unchanged. 
To further test the impact of the extinction value used 
to derive masses, we varied the Ay by a random amount 
in the interval ±2 mag from the spectroscopic value. We 
have done this 1000 times for each object, and finally 
adopted the average mas s fr om all the variations. The 
result is shown in Figure 10 with red open circles. We 
see that this has a small impact on the shape of the IMF. 
As mentioned earlier in the paper, members of Lupus 3 
typically have low extinction (< 5). In fact, the spec¬ 
troscopic sample used to plot the IMF contains 98% of 
members with Ay < 3. At Ay = 10, a O.O 7 M 0 object 
would fall right at the completeness limit of our survey, 
i.e. at the extinctions above 10 , we would not be able to 
detect substellar objects. However, we do not expect that 
many BD could be hiding at high extinctions, because in 
this case we should have also found some highly embed¬ 
ded stellar objects. The optical part of our spectroscopic 
survey would be sensitive to those, but we detect only 
one such source (SONYC-Lup3-l), whose spectrum and 
SED suggest a rare geometrical configuration (edge-on 
disk). 

The single power-law slopes of the mass functions for 
the two star forming regions derived here agree well, and 
are also consistent with typic al values found in the liter- 




atur' 

mass"stars and BD s below l ).6 Mf?^ in NGC 1333, whereas 


Scholz et al. (2012a) find a = 0.6 ± 0.1 for low- 


= 0.9 - I.tJ 


Oliveira et al. 
lor the substel 


Scholz et al. (2013), in a slightly modified analysis, report 
or th e same cluster. In IC348, |Alves d^ 


2013) report a of 1.0 ± 0.3 and 0.7 ± 0.4 
ar (< SOMjup) IMF extinction-limited. 


and spatially limited to the center of the cluster, respec¬ 
tively. This is c onsistent with a = 0 .7 — 0.8 found in the 


same c luster by Scholz et al. (2013). Alves de Oliveira 
et al.l (|2013|) also hnd a = 0.8 ± 0.4 and 0.7 ± 0.3 tor 


the substellar regime in p Oph, with the extinction lim¬ 
ited to Ay < 8 , and Ay < 15, respectively. Fo r low- 


mass stars and BDs in cr Ori (0.35 — O.OO 6 M 0 ), Pena 


Ramirez et al. (2012) find a = 0.6 ± 0 .2, in agreement 


power-law mass function, we first add the upper limit of 

2007 

Beiar et al. 12011 

). In Upper bco, assuming the 

11 missing members to the substellar bin. The objects 

ot 5-10 Myrs, Lodieu 

(2013 

); 

Lodieu et al. ( 

2013a), 


Ay=0 (0.009 M 0 ), and below the substellar limit. This 
case is extreme, both regarding t he n umber of missing 


port a = 0.5 — 0.7, tor t he ma ss range 


objects (the estimate from Section 4.2 was l'!l\°), and by 


the assumption that all the potentially missing members 
would be substellar, since some, or even all of them could 
equally be stellar. The result is shown in Figure rep¬ 
resented by blue crosses and histogram. After adding 
the missing members, the bins have been redistributed 


Zapatero Osorio et al. (2014) report a between 0 and 1 
lor the Pleiades cluster, down to « O.O 12 M 0 . 

We now outline the environmental conditions in our 
two clusters, to try to understand the origin of the dif¬ 
ference between the two IMFs in the substellar regime. 

Here we list only the most rec ently published results, for ear¬ 
lier works see [Scholz et al.'|||2012a[|. 
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TABLE 5 

Values for the mass functions show 
IN Figure!^ 


Mass interval 

Mean mass 

dN/dM 

Chamaeleon-I 

0.001 - 0.030 

0.0155 

465.52 

0.030 - 0.075 

0.0525 

344.44 

0.075 - 0.170 

0.1225 

168.42 

0.170 - 0.270 

0.2200 

130.00 

0.270 - 0.500 

0.3850 

58.70 

0.500 - 0.810 

0.6550 

40.32 

Lupus 3 

0.001 - 0.070 

0.0355 

210.15 

0.070 - 0.110 

0.0900 

375.00 

0.110 - 0.200 

0.1550 

188.89 

0.200 - 0.350 

0.2750 

83.33 

0.350 - 0.810 

0.5800 

23.91 


According to Gutermuth et al. (2009), with 19 stars pc“^ 
Lupus 3 has the lowest mean suriace density among 
nearby star forming regions, several times lower than 
e.g. pOph (120pc“^) and NGC 1333 (65pc“^), and 
also lower than Serpens (36pc“^), Gha-1 (30pc“^), CrA 
(29pc“^), or IC 348 (25pc“^). In terms of peak sur¬ 
face densities, Cha-1 is twice as d ense as Lupus 3 (800 vs 
400pc“^). In Scholz et al. (2013), we presented evidence 
for a difference m mass distributions of NGG 1333 and 
IG 348, under the (plausible) assumption that NGG 1333 
is closer to us. The denser of the two clusters (NGG 1333) 
harbors a larger fraction of BDs. Therefore, frequencies 
of BDs might depend on stellar densities, which is in line 
with predictions for BD formation through gravitational 
fragmentation of fila ments falling into a cluster potential 
( Bonnell et al.||2008[), dynamical cluster formation (sim- 
ulations by Bate 2012), or eve n the disk fragmen tation 
in some cases (see discussion in Scholz et al. ]|M^ . 

The presence of massive OB stars might influence BD 
formation in clusters, since BDs c ould form through 


photo-evapor ation of prestellar cores (Whitworth & Zin- 


necker 2004). However, nearby star forming regions 
(d 300 pc) contain extremely small populations of OB 
stars (only one B star in Lupus, and three in Cha-1), 
and their influence on the IMF of these clusters must be 
negligible. 


4.6. Brown dwarfs with masses below the 
deuterium-burning limit 

If we assume monotonic continuation of the power 
law shown in the upper panel of Figure 1^ below the 
D-burning limit, we would expect 8 ± 4 objects with 
0.005 < M < O.O15M0 in Cha-l. For the slope shown 
with the red dashed line, which accounts for the objects 
possibly missed by our follow-up, we get 12 ± 5. Our 
method for mass estimates from photometry yields 6 ob¬ 
jects with masses below ^ 0.015 Mq in the area of our 
survey, and we estimated that there might be 3 + 3 ° more 
in our candidate list. The data are therefore consistent 
with a monotonic power law mass spectrum across the 
D-burning limit, and down to ^O.OObM©. 

In NGC 1333, we find similar monotonic behavior of 
the low-mass IMF across the D-burning limit, but do not 
exclude a shallower slope below this limit (i.e. a < 0 . 6 ). 


In fact, based on more recent information about the mid- 
infrar ed colors of young ultracool objects (Faherty et al. 
2013), our upper limit on the number of the planetary- 


mass BDs is probably too conservative, and is likely to 
be lower (4 instead of 8 ). This suggests a possible break 
in the power-law arou nd D-burning limit. In a Ori, |Pefia] 
Ramirez et al. (2012|) report the IMF consistent with a 
smooth power law with a = 0.6 ± 0.2 down to O.OOdM©. 
In the range 0.004 - O.OO 3 M 0 there are fewer candidate 
members observed in their deep J-band catalog than one 
would expect from an extrapolation of the same power- 
law to these mass es, possibly signaling a lower slope in 
this mass regime. Lodieu et al. (2013a) present a deep 
NIR survey probing the planetary-mass regime in UpSco, 
complementin g their previous work in the same region 
(|Lodieu||2013l). The updated IMF in UpSco, extending 
down to ~ O.OO 5 M 0 , is consistent with a rising power-law 
function with a ~ 0.5. 

In the area of our Cha-1 survey, there are in total 
94 known members, out of which ~6 appear to have 
masses below O.O15M0. Previously we estimated that 
there might be 0-13 objects missing in this mass regime, 
and 4-13 among the higher mass BDs. In the current cen¬ 
sus, the brown dwarfs with masses in the planetary-mass 
regime seem to comprise only ^ 6 % of the Cha-1 pop¬ 
ulation, whereas their contribution can go up to ^17% 
taking into account these estimates. Even in the extreme 
case, the mass budget of these objects can be only about 
1% of the total mass of the cluster. 

In Lupus 3, our survey is sensitive only to the high- 
mass tip of the planetary-mass regime, but |Comerdn| 
( 2011 ) conducted a small-area photometric survey sev- 
erai rnagnitudes deeper than ours and sensitive to free- 
floating Jupiters according to the models. The saturation 
of their survey matches the completeness limits of ours, 
and their 3a detection limit is at Iq = 25.6 mag. The 
number of the cool photometr ic candidates identified in 
the survey by Comeron (2011) is consistent with statis¬ 
tical expectation of a background population. The sur¬ 
vey encompasses a very small part of Lupus 3, about 100 
times smaller than the SONYC area, but carefully chosen 
to match an area known to be rich in higher mass stellar 
and substellar members. While a wider area survey sen¬ 
sitive to t he lowest masses is still required to confirm the 
result by Comeron (2011), it represents a complement 
to our finding that there are very few substellar mem¬ 
bers left to be discovered down to the mass limits of our 
survey. 

To conclude, surveys in various star forming regions 
clearly show that the IMF below ~O. 6 M 0 , and down to 
^O.O 2 M 0 can be well described by a monotonic power- 
law with an exponent a = 0.6—1.0. Below the D-burning 
limit, the slope seems to be similar, or shallower, but it 
is certainly not steeper than the one ab ove 0 . 02 Mf 7 >. This 
means that, if the microlensing results (Sumi et al.|2011 ) 
hold, the objects they probe must undergo a different 
formation path than those studied so far in young clus¬ 
ters. 


5. CONCLUSIONS AND SUMMARY 

Here we have presented the further spectroscopic 
follow-up of the VLM candidate member lists in Cha- 
1 and Lupus 3, using the NIR spectrograph SofI at the 
NTT. This work is a continuation of our previous efforts 
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described in Muzic et al. (2011 1 and Muzic et al. (2014). 

In Lupus 3, we obtained 19 new spectra ol our high pri- 
ority photometric and proper motion sample (“iJ-pm”), 
thus almost completing the spectroscopic survey of this 
sample: we have targeted 50 out of 53 candidates above 
the surveys completeness limit. We identified two prob¬ 
able substellar members of Lupus 3, one of which is new 
and has a NIR spectral type M7.5. Based on statistical 
arguments, we conclude that only a very few substellar 
members are left to be discovered in this region, at least 
down to the completeness limit of our survey, lying at 
i = 20.3, equivalent 0.008 — 0.02M 0 , for Ay <5 at 1 
Myr. 

The follow-up in Cha-I included 15 candidates with 
expected masses in the planetary regime, of which only 
one was confirmed as substellar with the spectral type 
L3 and Teg of 2200 K. According to the BT-Settl mod¬ 
els, an object at this Teg should have a mass of ^ 
0.009 — O.OI 2 M 0 at the age of 1-2 Myr. The compar¬ 
ison of the multi-wavelength photometry of this object 
with the models yields a mass of ^ 0.007 M 0 , at the dis¬ 
tance of Cha-I and age of 2 Myr. Based on statistical 
arguments, we estimate that in the area of our survey, 
Sllil objects are left to be discovered with 17.5 < / < 21, 
and 3 ^ 3 °% with 21 < / < 23. In terms of masses, the 
first bin is equivalent 0.005 — 0.06 M 0 , and the second 
one 0.003 — 0.02 M 0 , at the age of 1-2 Myr, and Ay < 5. 


The IMF below IM 0 in both Cha-I and Lupus 3 can be 
described by a power law dN/dM oc M““, with the slope 
a ^ 0.7, in agreement with the results of the surveys in 
other clusters. The same is valid for the star-to-BD ra¬ 
tio, which is found to be between 1.5 and 6 . In Lupus 3, 
however, we find evidence for a flattening of the slope, 
or even a possible turnover of the IMF in the power law 
form in the substellar regime: this region seems to pro¬ 
duce less BDs in comparison to majority of other clusters. 
The flattening is still present even after accounting for 
the maximum number of objects possibly missed by our 
survey. 

The IMF in Cha-1 shows a monotonic behaviour across 
the D-burning limit. The expected total number of sub¬ 
stellar objects can be statistically estimated from the 
success rates of our spectroscopic follow-up, and this es¬ 
timate yields numbers consistent with the same power 
law extending down to our completeness limit, lying at 
0.004 - 0.009 M 0 . We estimate that the low-mass mem¬ 
bers below the D-burning limit contribute of the order 
5 — 15% to the total number of Cha-I members, compris¬ 
ing therefore < 1 % of the mass budget of the cluster. 


The authors thank Katelyn Alters, Catarina Alves de 
Oliveira, and David Lafreniere for sharing their spectra, 
used for comparison in this work. This work was co¬ 
funded by NSERC grants to RJ. 
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